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LONG-LIVED MAGNETIC-TENSION-DRIVEN MODES IN A MOLECULAR CLOUD 
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ABSTRACT 

We calculate and analyze the longevity of magnetohydrodynamic (MHD) wave modes that occur in 
the plane of a magnetic thin sheet. Initial turbulent conditions applied to a magnetically subcritical 
cloud are shown to lead to relatively rapid energy decay if ambipolar diffusion is introduced at a level 
corresponding to partial ionization primarily by cosmic rays. However, in the flux-freezing limit, as 
may be applicable to photoionized molecular cloud envelopes, the turbulence persists at "nonlinear" 
levels in comparison with the isothermal sound speed c s , with one-dimensional rms material motions 
in the range of ~ 2 c s — 5 c s for cloud sizes in the range of ~ 2 pc — 16 pc. These fluctuations 
persist indefinitely, maintaining a significant portion of the initial turbulent kinetic energy. We find 
the analytic explanation for these persistent fluctuations. They are magnetic-tension-driven modes 
associated with the interaction of the sheet with the external magnetic field. The phase speed of such 
modes is quite large, allowing residual motions to persist without dissipation in the flux- freezing limit, 
even as they are nonlinear with respect to the sound speed. We speculate that long-lived large-scale 
MHD modes such as these may provide the key to understanding observed supersonic motions in 
molecular clouds. 

Subject headings: instabilities — ISM: clouds — ISM: magnetic fields — magnetohydrodynamics 
(MHD) — stars: formation 



1. INTRODUCTION 

Nont hermal linewidths a re ubiquitous in molecular 
clouds (|Solomon et al.|[l987j) and are interpreted to rep- 
resent highly supersonic random internal motions (see 
iMcKee fc Ostrikeri[2007l for a recent r eview) . Principal 
component analysis (jBrunt et al.ll2009t ) reveals that most 
of the energy is contained in modes that span the largest 
scale of the cloud. 

Since molecular clouds are threaded by large-scale 
magnetic fields, an attractive suggestion has been that 
the turbulence represents supersonic but sub-Alfvenic 
magnetohydrodynamic (MHD) waves, with the noncom- 
pressive shear Alfv en mode ident i fied a s a possible long- 
lived component ()Arons fc Maxl [19751) . This was in- 
tended to bypass the usual problem of rapid dissipation 
of supersonic hydrodynamic turbulence through shocks. 
A compilation of available Zeeman measurements of the 
(line-of-sight) magnetic field strength, gas density p, and 
one-dimensional velocity dispersion a v shows that cloud 
fragments obey a direct linea r correlation between <j v 
and the me an Alfven speed Va (|Mvers fc Goodmanlll988t 
iBasul 120001 ), with a proportionality coefficient ~ 0.5. 
This correlation may be attributed to a rough equality of 
the magnitudes of gravitational, magnetic, and turbulent 
energi es, a nd was interpreted by IMouschovias fc Psaltisl 
(| 19951) and IMouschovias et al.l (|2006f ) to mean that the 
motions are Alfvenic disturbances in which the perturbed 
magnetic field is comparable in strength to the back- 
ground magnetic field. Furthermore, the MHD mo- 
tions in a molecular cloud ma y represent long- wavelengt h 
standing waves, as argued bv IMouschovias! (|1975l Il987f l. 
This brings out the possibility of "global" effects (e.g. 
due to cloud boundaries and external interaction) being 
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important in understanding cloud turbulence, and the 
need to go beyond comparing observations with models 
of wave propagation in an infinite medium. 

In direct contrast to the scenario of long-lived mo- 
tions, numerical simulations of molecular cloud turbu- 
lence using a three-dimensional simulation cube with pe- 
riodic boundary conditions have revealed that supersonic 
MHD turbulence decays away rapidly, like its hydrody- 
namic counterpa rt, on about a sound crossing time of 
the driving scal e dStone et al.lll998l:lMac Low et al.lll998t 
iMac Low! [19991 : lOstriker et all 120011 ). This happens in 
either the case of sub-Alfvenic or super-Alfvenic turbu- 
lence, and in both cases, turbulence is maintained for 
long periods only by constant driving of velocity pertur- 
bations in Fourier space. When interpreting the above 
results, we should keep in mind that periodic box simu- 
lations represent a "local" patch of uniform background 
density that is embedded within a larger cloud, and are 
equivalent to studying an infinite uniform medium. By 
comparison, a 1.5-di mensional global model inc luding 
vertical stratification (jKudoh fc Basuir2003l [2006) found 
that the decay of turbulence could be delayed, but only 
mildly, by some transfer of internal kinetic energy from 
small to large scale modes along the magnetic field di- 
rection. The rapid turbulence dissipation in all of these 
models is due to the presence of shocks and takes place 
under the assumption of magnetic flux freezing, without 
any contribution from magnetic field dissipation, e.g., by 
ambipolar diffusion. 

The bottom line from the above studies is that all 
previous numerical modeling of MHD turbulence leads 
to rapid dissipation, in about a crossing time, with a 
logical conclusion that matching observations requires 
vigorous driving of turbulence from unspecified sources. 
The alternate possibility of maintaining some long-lived 
global modes is appealing but remained largely unex- 
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plored quantitatively. 

In a recent paper, iBasu et all (|2009bl ) carried out 
an extensive parameter survey of fragmentation initi- 
ated by nonlinear turbulent flows, employing the mag- 
netic thin-sheet approximation and als o including the 
effect of ambipolar dif f usion ( see also IBasu fc Ciolekl 
20041: ILTfc Nakamural [200l jNakamura fc Lil 120051: 
Ciolek fe Basull2006t IBasu et al.H2009al) . In this approx- 
imation, the sheet interacts at its upper and lower sur- 
faces with an external magnetic field, and can be con- 
sidered a global model in the z-direction (parallel to the 
mean background magnetic field), although it is a local 
period ic) model in the x- and y-directions. IBasu et al.l 
2009b|) found that initial turbulent fluctuations decayed 
away quite rapidly in all models with supercritical mass- 
to-flux ratio, as well as for subcritical models that in- 
cluded the effect of ambipolar diffusion. However, a sur- 
prising result was that subcritical clouds evolving under 
flux-freezing were able to maintain a substantial portion 
of their initial kinetic energy to indefinitely large times. 

In this paper, we analyze this unique instance of a tur- 
bulent MHD simulation that yields long-lived nonlinear 
motions. We perform a suite of numerical simulations to 
test its generality, and also establish an analytic expla- 
nation for this very interesting result . 

2. METHOD 

The thin-sheet equations are obtained by vertically in- 
tegrating the full equations for a partially ionized, mag- 
netic, self-gravitating, isothermal fluid in the vertical di- 
rection from z = — Z{x,y) to z = +Z(x,y). Details 
of thi s inte gration are found in ICiolek fc Mouschoviasl 
(|1993fl and ICiolek fc Basul (|2006f ). The nonaxisymmet- 
ric thin-sheet equations, formulation of our model, and 
our numerical methods are described in several papers 
(ICiolek fc Basul 120061: IBasu et al.l [2009al lbT). The evo- 
lution equations for the magnetized thin sheet include 
the effects of magnetic tension due to the external mag- 
netic field B(x 7 y, z). It is calculated as a potential field, 
with the vertical magnetic field in the equatorial plane, 
B z ,eq{x, y), acting as a source for B(x, y, z) — B re [Z, much 
as a n (x, y) acts as a source for the gravitational field. Pe- 
riodic boundary conditions are applied to a square Carte- 
sian region of size L. The initial background state has a 
uniform (in the x- and y- directions) neutral surface den- 
sity i7n,o and a uniform vertical (z-direction) magnetic 
field B Te f. We input nonlinear velocity fluctuations with 
spectrum v\ oc k~ A in Fourier space, where k is the abso- 
lute value of the wavenumber, and modes are damped at 
a fixed (small) scale that is independent of the box size 
or the number of grid zones used in a simulation. 

The gas is isothermal with sound speed c s , and par- 
tial ionization is mainly due to cosmic rays. This intro- 
duces the dimensionless free parameter fn^o = TV^n/to, 
where r^o is the initial neutral-ion collision time, and 
to = c s /2i:G<T n fl is a characteristic time in the prob- 
lem. The flux-freezing limit, used extensively in this 
paper, corresponds to f^o = 0. Another important 
parameter is the initial dimensionless mass-to-flux ratio 
Ha = 2nG 1 ^ 2 a n fi/ B vc {, i.e., /io > 1 yields a supercriti- 
cal cloud in which fragmentation occurs dynamically and 
Ho < I yields a subcritical cloud in which fragmentation 
is d riven on a longer tim e scale by ambipolar diffusion 
(see ICiolek fc Basul 120061) . Turbulent initial conditions 
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Fig. 1. — Kinetic energy Bkin normalized to its initial value Ekin.o 
for two models. Each model has no = 0.5, v a = 2c s , and L/Lo = 
327r. One model evolves with flux freezing, f a i q = 0, and is run 
with N = 512. The other model has partial coupling of neutrals 
and ions characterized by f n j o = 0-2, and is run with N = 1024. 

introduce the dimensionless parameter v a /c s , where v a 
is the rms amplitude of the initial velocity fluctuations 
in each direction. Finally, we also vary the ratio L/Lq 
(where Lq = c 2 /27rGc7 ni o is a characteristic length scale 
of the system), and N, the number of grid zones in each 
direction. An additional parameter, the dimensionless 
external pressure P cx t = 2P cx t/7rGcr^ , is kept fixed at 
0.1 in all models, and does not play an important role in 
the dynamics. 

3. RESULTS 
3.1. Canonical models 

Two models illustrate the key result for subcritical 
clouds with turbulent initial conditions. Figure [T] shows 
the time evolution of total kinetic energy for a model 
which allows for neutral-ion drift (t^o = 0.2), and an- 
other model which has flux freezing (f n i,o — 0). The 
value fn^o = 0.2 corresponds to the canonical ioniza- 
tion fraction implied by primarily cosmic ray ionization: 
Xi ^ 10~ 7 (n n /10 4 cm" 3 )- 1 / 2 (see iTielensI [2001 . where 
n n is the number density of neutrals. Both models are 
characterized by fj,o = 0.5, v a = 2c s ,L = 32ttLq. The 
flux- frozen run has N = 512 while the ambipolar diffu- 
sion run has N = 1024. The evolution of the ambipolar 
diffusion model terminates at time t = 45.4 to, when the 
highest column density in the simulation reaches 100 a n> o 
— a useful proxy for runaway collapse of the first core. 
At this time, the kinetic energy has decayed away sub- 
stantially, and appears to still be declining. In contrast, 
the flux- frozen model has, after an initial loss of some 
kinetic energy, stabilized to executing oscillations about 
a mean value -Ekin ~ 0.7£'kin,o- We have run simulations 
with flux freezing up to t ss 35, 000 to with no change in 
this behavior. 

Figure [2] shows color images of the column density for 
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Fig. 2. — Images of gas column density cr n (x, y)/(?n,o for models 
with ambipolar diffusion (top) and flux freezing (bottom), shown 
in identical color schemes that are proportional to the logarithm of 
the column density. Both models have [iq = 0.5, initial turbulence 
with v a = 2 c s and spectrum u? oc k~ 4 , and are run with N = 256. 
The ambipolar diffusion model terminates at time t = 43.5 to in 
this realization due to the eventual runaway collapse of a core (in 
the upper right of the image). The flux-frozen model is shown at 
the same time, but continues to evolve to indefinitely large times 
without collapse. It shows a more wispy column density structure, 
with no evidence of monolithic collapse toward any density peaks. 
An animation of the evolution of each model is available online. 

models that are equivalent to the ones described above, 
but with N = 256. The ambipolar diffusion model is 
shown at its end time t = 43.5 t Q , and the flux-frozen 
model is shown at the same time, although it contin- 
ues to evolve indefinitely. The ambipolar diffusion model 
shows evidence of monolithic collapse toward one or more 
density peaks, while the flux-frozen model shows a more 
wispy character and gives no indication of impending col- 
lapse in any region, neither visually nor quantitatively. 
Figure |3] shows a model snapshot of the flux-frozen model 
but viewed from a three-dimensional perspective, with 
the external field lines illustrated in the region above the 
sheet. Since this is a subcritical model, the field lines are 
not significantly deformed. The pitch angle of the mag- 
netic field relative to the vertical direction, measured at 
the sheet surface, maintains an average value that is a lit- 
tle less than 10°. For a compariso n of field line curvat ure 
for models with a range of /x , see lBasu et al.l (|2009al lbh. 
Animations of both Figures [2] and G3 with the latter also 
showing external field line evolution, are available online. 




Fig. 3. — Surface map of gas column density a n (x, y)/a n ,o for 
the same flux-frozen model that is shown in Fig. [2] The eleva- 
tion and color of the surface is proportional to the logarithm of 
the local column density. The sheet is viewed from a non-face-on 
viewing angle. The external magnetic field above the sheet is also 
represented by field lines, in a region near a density peak. An ani- 
mation of the evolution of the model, including field line evolution, 
is available online. 

3.2. Connection to Linear Analysis 

To gain insight into the long-lived mode for the flux- 
freezing model, we revisit the modal analysis of a par- 
tially i onized magnet ized sheet. This ha s been carried 
out bv lMortonl (|1991f ) and lCiolek fc Basul (|2006[ ). Start- 
ing with a static uniform background, the thin-sheet 
equations can be expanded by writing any physical quan- 
tity f(x,y,t) = fo + S f a exp[i(k x x + k y y - ut)], where 
/o is the unperturbed state, Sf a is the amplitude of the 
perturbation, k x and k y are the wavenumbers in the x- 
and y- directions, and u) is a complex angular frequency. 
For assumed small perturbations, one can retain only 
terms up to first order in pertu rbed quantities, result- 
ing in Equations 32(a) - 32(d) of iCiolek fc Basul ([2006). 
Those equations can be combined to yield the dispersion 
relation 

(u + id) [u 2 - C 2 s k 2 + 2nGa nfl k] 

= Lu[Vl Q k 2 + 2TrGa n Wv 2 k] , (1) 

where 

e = 2vr r ni , [ Vl fl k 2 + 2nGa nfi nfk] (2) 

contains the effect of ambipolar diffusion. In the above 
equations, we have introduced the Alfven speed Va,o for 
physical clarity of the magnetic-pressure-driven terms 
proportional to k 2 (= k 2 + k 2 ), while retaining /j,q for 
physical clarity of the magnetic-tension-driven terms pro- 
portional to fc (= |vA?|)- The two parameters are ac- 
tually related: F| = B 2 ci /4:Trp nfi = 2nGa nfi ^ '' Z o, 
where the sheet half-thickness Zq = a n . o/2p n ,o and p n< o 
is its mass volume density. The quantity C c ff is an ef- 
fective sound speed that takes into accou nt the restoring 
force due to an external pressure P ex t (see lCiolek fc Basul 
2006). In the flux- freezing limit, Equation ([T]) becomes 

^ = (Vl, a + C 2 S ) k 2 + 27rGcr n , (p Q 2 - l)fc. (3) 

In the case of subcritical clouds (fi < 1), the second 
term on the right-hand side becomes a stabilizing term 
rather than a destabilizing term associated with gravi- 
tational instability. However, the full dispersion relation 
(Equation ([T])) does contain destabilizing terms due to 
ambipolar diffusion; the effect on a subcritical cloud is 
a "slow" instability leading to collapse on an ambipo- 
lar diffusion tim e scale rather t han a dynamical time. 
Equation (42) of lCiolek fc Basul (|2006l) yields a good ap- 
proximation to the ambipolar diffusion growth time for 
significantly subcritical clouds. 
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Equation © shows that long- wavelength modes evolv- 
ing under flux freezing have a phase speed 



MT,0 



(/V 



(4) 



where A = 2n/k. These modes are driven by the restor- 
ing force of the magnetic tension of inclined field lines 
that connect the sheet to the external medium. These 
magnetic-tension-driven modes should not be confused 
with the traditional MHD wave modes. Within the 
thin sheet, in the short-wavelength limit, magnetic pres- 
sure drives the magnetosound mode with phase speed 

^MS,0 = Oto + 1/2 - 

Since Vmt,o vA, it achieves significant values (much 
larger than Vms,o)j f° r Mo = 0-5 and wavelengths equal 
to the box sizes we consider: L = 16tt Lq — 1287rLo- The 
values are in the range 4.9 c s — 13.9 c s and typical values of 
input parameters wou ld correspond to d imensional box 
sizes « 2 — 16 pc (see lBasu et al.ll2009al lb1. for scaling of 
units). Since the restoring force is provided by the ex- 
ternal potential field that can adjust instantaneously as 
the sheet evolves, the modes found in this linear analysis 
cannot be applied to arbitrarily large wavelengths. In 
reality, there must be time for readjustment of the exter- 
nal field. The magnetic potential ^ m(x, v, z) above and 
below the sheet decays a s exp (— k\z\) (see lCiolek fc Basul 
120061 : IBasu et al"]|2009a|) . so that a characteristic height 
of deformation of the field lines is fc _1 . The ratio e of 
the Alfven crossing time across this distance divided by 
the wave period must be well below unity in order for 
the potential field approximation to be valid. While the 
Alfven crossing time grows more rapidly (oc A) with in- 
creasing wavelength than does the wave period (oc \/A), 
we find that e < 0.26 for modes of even our largest box 
size, if the external density p cx t < 0.1 p n ,0- The nature 
of a low-density medium external to clouds or clumps is 
discussed in Section [4] 

An interesting analogy can be made between the 
magnetic-tension-driven modes and gravity-driven waves 
in deep water. There, the undulations of wavenumber k 
on a water surface can be felt down to a characteristic 
depth k~ x . Velocities below the surface are determined 
from a velocity potential solution of Laplace's equation. 
This is partly due to the incompressible fluid approxi- 
mation in which the water pressure can adjust instan- 
taneously. A clearer mathematical analogy also occurs 
in the following manner. Since the vertical gravitational 
field above and below a uniform thin sheet has magnitude 
\g z \ = 2irGa n fi, Equation (j4]) may be rewritten as 



Vmt,o = vVo 2 - > 



(5) 



in analogy to the phase speed v — \J gjk for deep water 
waves in a constant gravitational field g. 

3.3. Further modeling 

The implication of the high values of Vmt.o for the 
magnetic-tension-driven modes in the plane of a thin 
sheet is that waves with nonlinear particle motions (in 
comparison to the sound speed c s or magnetosound speed 
Vms,o) may still act as linear waves since their material 
motions are much slower than Vmt,o- They will then 



evolve (in the flux-freezing limit) without any nonlinear 
distortion and dissipation. All models do lose signifi- 
cant kinetic energy in an early phase, due to shocks and 
compression that leads to significant losses in an isother- 
mal gas. Small-scale modes are prone to such decay; for 
them the relevant signal speed is Vms,o> which equals 
2.9 c s for the models presented in Section IBTTl However, 



the 



k spectrum means that most of the energy 



is in the largest scale mode, which can survive indefi- 
nitely, as long as the velocity amplitude is significantly 
less than Vmt,0j either initially or after some nonlinear 
decay. In the magnetic-tension-driven mode, energy is 
stored and released by the magnetic field, without losses 
due to ambipolar diffusion (in the f^o = models) or 
other forms of magnetic field dissipation. Furthermore, 
the isothermal assumption does not rob any net energy 
at this stage. In symmetric oscillations, energy is lost 
during wave compressions and an equal amount gained 
back during wave expansions. 

Figure [4] explores the effect of different initial condi- 
tions on the decay and residual amount of turbulence 
in several flux-frozen (t^o = 0) models. The top 
panel shows a comparison of models with vf. proportional 
to fc~ 4 , fc~ 2 , and k°, respectively, but all having the 
same initial rms speed. The spectrum with the greatest 
amount of energy on the largest scale retains the most en- 
ergy, as it is the largest scale mode that has the greatest 
phase speed and is most likely to survive with significant 
amplitude. The bottom panel shows that models with 
fixed spectrum v\ oc fc -4 but differing v a will lose dif- 
ferent proportions of their initial turbulent energy. The 
phase speed of the largest mode in these simulation boxes 
is Vmt,o = 4.9 Cg, and increasing v a leads to greater pro- 
portionate loss of initial kinetic energy. However, there is 
a weak trend toward retaining a greater absolute amount 
of energy, as a tabulation of v xma , the one-dimensional 
rms speed at t ~ 100 to in each simulation, reveals. Ta- 
ble [1] lists v rms for many models that have flux- freezing, 
/io = 0.5, and v\ oc fc -4 initially. Supersonic motions re- 
main in all models, and the residual amplitude rises with 
increasing box size as well as initial velocity amplitude 
v a . The values of w rms appear to saturate however, so 
that they remain a reasonably small fraction (14%-56%) 
of Vmt.o for each model. 

Oscillations of the kinetic energy are clearly visible in 
the models that retain a large part of their initial en- 
ergy, so that the values of v ims in Table Q] are varying by 
up to 10%. The dominance of the largest mode in the 
initial conditions and the preferential nonlinear damping 
of smaller modes implies that the period of the largest 
mode is a reasonable approximation to these observed 
periods P. We determine P by an average over many 
peak-to-peak oscillations in each model. The kinetic en- 
ergy should oscillate with half the period of the velocity 
cigcnfunction, so its expected period is 



P = -- 



1 L 



2 %T,0 



L 



(fi - l)Gcr n ,o 



1/2 



(6) 



where we have used A = L. In terms of the dimensionless 
box size L/Lq and to, we can write this as 



1 

P =2 



2n(L/L ) 



Mo 2 - 1 



1/2 



to- 



(7) 
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Fig. 4. — Effect of different power spectra and turbuient velocity 
amplitudes on kinetic energy decay. All models evolve with flux 
freezing (-f n ; q = 0). Top panel: three models with fixed v a and 
other parameters (as labeled) but differing power spectra of initial 
turbulence, as labeled next to each curve. The model with the most 
energy on the largest scale retains the greatest part of its energy. 
Bottom panel: three models with differing v a but all other param- 
eters including the power spectrum kept fixed, as labeled. Models 
with greater v a lose a greater proportion of their initial kinetic 
energy, although they retain a slightly greater absolute amount of 
kinetic energy (see Table [TJ. 



TABLE 1 

Velocity Amplitude Results for Selected Models 



N 


L/L 


Vmt,o/ c s 


Va/c B 


"rras / C-s 


l>rms/VMT,0 


128 


16?r 


4.9 


2 


1.7 


0.35 


128 


16tt 


4.9 


i 


2.4 


0.48 


128 


16tt 


4.9 


6 


2.6 


0.53 


256 


32tt 


6.9 


2 


1.6 


0.23 


256 


32tt 


6.9 


1 


3.4 


0.49 


256 


32tt 


6.9 


6 


3.9 


0.56 


512 


32tt 


6.9 


2 


1.7 


0.25 


512 


64tt 


9.8 


2 


1.9 


0.19 


512 


64tt 


9.8 


3 


2.7 


0.28 


512 


64tt 


9.8 


1 


2.9 


0.29 


1024 


128tt 


13.9 


2 


1.9 


0.14 


1024 


1287T 


13.9 


4 


3.4 


0.25 


1024 


1287T 


13.9 


6 


5.0 


0.36 



Note. — All models above have fj,o = 0.5, f n ifi = 0, and initial 
turbulent spectrum oc k~ 4 . The one-dimensional velocity dis- 
persion » rms is measured at t 100 to and is present indefinitely 
with variability of up to 10%. 



Figure [5] shows the predicted dependence in solid lines, 
for fi — (0.25,0.5,0.7). Different symbols as described 
in the figure caption represent the empirical determina- 
tions of P from various models. The agreement is re- 
markably good, and improves for the largest box sizes, 
where the long- wavelength approximation made in Equa- 
tion (01 holds particularly well. 

4. DISCUSSION 

The study of the decay of MHD turbulence has gener- 
ally been based on the modeling of Alfven, slow MHD, 
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Fig. 5. — Comparison of analytically predicted periods of kinetic 
energy oscillation with results of simulations. Solid lines are the 
predicted periods from the linear theory of large-scale flux-frozen 
modes driven by magnetic tension, for three different values of 
the dimensionless mass-to-flux ratio fio = (0.25,0.5,0.7). Open 
diamonds represent average periods of oscillation for four models 
of differing box size L and fixed fio = 0.7. Asterisks and plus 
signs represent the same but for fixed fio = 0.5 and fio = 0.25, 
respectively. The agreement is between the solid lines and symbols 
is remarkably good, and best for the largest box sizes, since the 
analytic prediction is made in the long- wavelength limit. 

and fast MHD modes in media that have a uniform back- 
ground. The study of more complex global MHD (includ- 
ing magnetogravitational) modes for molecular clouds re- 
mains to be explored. In this paper, we have analyzed 
turbulent decay in a magnetically subcritical sheet-like 
cloud. It is an idealized model of a molecular cloud that 
is tied to a magnetic field anchored in the interstellar 
medium. A large fraction of the initial input kinetic en- 
ergy is retained by the deformed magnetic field, and then 
persists in the cloud as large-scale oscillations. These 
represent linear waves of large extent which are never- 
theless supersonic since the phase speed of the magnetic- 
tension-driven modes is up to ~ 10 times the sound speed 
for typical cloud sizes. 

Our model may approximate the situation of molec- 
ular clouds that are embedded in a low-density warm 
H I halo, or even the case of molecular cloud clumps 
that may be embedded in a matrix of H I gas (see 
iHennebelle fe I nutsuka 200(f) . Flux freezing is a good ap- 
proximation for molecular cloud envelopes (as opposed to 
dense cores), due t o significant photoionization by back- 
groun d starlight (jMcKeel 119891: iCiolek fe Mouschoviasl 
1995). Observations also reveal that the low-column- 
density molecular clo ud envelopes actually contain most 
of the cloud mass (jKirk et all I200& iGoldsmith et~aTl 
2008). These envelopes may have a subcritical mass- 
to-flu x ratio, as impl ied by their lack of star forma- 
tion (iKirk et al.l 2006D . velocity data (e.g., in Taurus, 



iHever et all I200S), and the sub critical state of the H 
I clouds (jHeiles fe Trolandl l2005h from which molecular 
clouds are presumably assembled. 
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Continuous driving of turbulent motions in molecular 
clouds is often invoked because the canon ical numeri- 
cal result of decay in a crossing time (e.g., iStone et"al] 
119981 : iMac Low! 119991 ) is inconsistent with estimated 
cloud lifetimes that are at least a few cr o ssing times 
(| Williams fc McKed [19971) . IBasu fc Murahl (pOOll ) have 
argued that continuous driving of turbulence is consis- 
tent with observational constraints only if the driving 
occurs on the largest scale in the cloud, i.e., most of 
the energy is contained on that scale. Furthermore, con- 
tinuous driving may not even be required if the decay 
time of the large-scale modes is greater than or equal to 
the estimated cloud lifetimes. Our models suggest that 
large-scale modes that are coupled to the external mag- 
netic field can persist for very long times, thus reducing 
the need for continuous driving in order to explain ob- 
servations. These modes preferentially span the largest 
scales in the model cloud, i n agreement with analysis 
of observed cloud turbulence (jBrunt et al.ll2009l) . Future 
spectral line modeling of the large-scale cloud oscillations 
in our model cloud may make for useful comparison with 
observations, as has been done in a previous study o f 
motions in the vicinity of dense cores (|Kirk et alJ [2009). 

A counter-effect to the maintenance of large-scale 
modes is the loss of energy to the external medium. 
This can be accomplished by a coupling of the magnetic- 
tension-driven modes to MHD modes in t he external 
medium. In a related example, lEnd (|2002f ) found sig- 
nificant energy loss to an external medium during core 
contraction using an approximate treatment of the trans- 
mission of transverse Alfven waves through the bounding 
surfaces of a thin sheet. Some form of MHD wave cou- 
pling is certainly at work between any molecular cloud 
and its environment, although one may also argue that a 
clump embedded in a larger complex may reach a steady 
state in which it gains as much energy from its exterior 
as it loses. In any case, the study of the propagation 
of waves outside the cloud is outside the scope of our 
model. Future three-dimensional global MHD models of 



molecular clouds, which include the effect of an external 
medium, can address this point. 

5. SUMMARY 

We have demonstrated that MHD modes driven by the 
tension of inclined magnetic held lines have a large phase 
speed for subcritical clouds, which increases in propor- 
tion to the square root of the wavelength. Numerical sim- 
ulations show that nonlinear motions (in comparison to 
the sound speed) persist indefinitely for thin-sheet evolu- 
tion in the limit of magnetic flux-freezing. These are the 
first of any variety of MHD turbulence simulations that 
show long-lived nonlinear motions. For the broad set of 
models that we study, the residual one-dimensional rms 
material motions are in the range of w 2 c s — 5 c s for cloud 
sizes in the range of « 2 pc — 16 pc. We find that run- 
away collapse toward isolated density peaks occurs when 
partial ionization due to (primarily) cosmic rays and am- 
bipolar diffusion is included. However, the flux-freezing 
results can be relevant to understanding the low-column- 
density molecular cloud envelopes, which are photoion- 
ized to the level of effective flux-freezing, and contain 
most of the mass in a molecular cloud. For those regions, 
there is a particularly important role for wave modes 
driven by a magnetic field that threads the cloud and 
is connected to an external medium. Long-wavelength 
modes such as the ones we study may provide at least 
part of the explanation for widely observed supersonic 
motions in molecular clouds. 
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